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I ■ ABSTRACT 

O^l \ We present the first detection of a gravitational depletion signal at near-infrared wave- 

lengths, based on deep panoramic images of the cluster Abell 2219 (z=0.22) taken with 
the Cambridge Infrared Sm'vey Instrument (CIRSI) at the prime focus of the 4.2m 
William Herschel Telescope. Infrared studies of gravitational depletion offer a number 
of advantages over similar techniques applied at optical wavelengths, and can provide 
reliable total masses for intermediate redshift clusters. Using the maximum likelihood 
. technique developed by Schneider, King & Erben (1999), we detect the gravitational 

' depletion at the 3ct confidence level. By modeling the mass distribution as a singular 

\ isothermal sphere and ignoring the uncertainty in the unlensed number counts, we find 

i an Einstein radius of 0e — 13.7^42 arcsec (66% confidence limit). This corresponds 

\ to a projected velocity dispersion of ct^, ~ 800 km s~^, in agreement with constraints 

pi • from strongly-lensed features. For a Navarro, Frenk and White mass model, the ra- 

I ' dial dependence observed indicates a best-fitting halo scale length of 125 h~^ kpc. We 

O \ investigate the uncertainties arising from the observed fluctuations in the unlensed 

-4-^ ■ number counts, and show that clustering is the dominant source of error. We extend 

^ ' the maximum likelihood method to include the effect of incompleteness, and discuss 

. . , the prospects of further systematic studies of lensing in the near-infrared band. 
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1 INTRODUCTION 

Gravitational lensing is now recognised as a valuable probe 
of the mass distribution in intermediate redshift {z < 0.4) 
galaxy clusters independent of any assumptions about the 
nature of the lensing material (see reviews by Fort & Mel- 
lier 1995, Narayan & Bartelmann 1997, Mellier 1999, Bartel- 
mann & Schneider 1999). For lenses at such redshifts, uncer- 
tainties arising from the redshift distribution of background 
sources are minimized and the angular scales of both weak 
and strongly-lensed features are well-suited for precise stud- 
ies. 

On small angular scales in super-critical systems, 
multiply-lensed arcs can provide useful absolute mass es- 
timates provided spectroscopic redshifts and geometric con- 
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straints on the location of the various critical lines are avail- 
able (e.g. Kneib et al. 1994, 1996). The most promising 
progress in constraining the mass on large scales has come 
from weak shear measurements (e.g. Fischer 1999; Hoekstra 
et al. 1998; Clowe et al. 1998) for which sophisticated inver- 
sion techniques have been developed (e.g. Kaiser & Squires 
1993; Kaiser 1995). 

However, cluster mass determinations based on weak 
shear signals are not without limitations. Firstly, the mass 
reconstruction from weak shear is non-trivial because of 
boundary effects due to the finite field of the data. To 
counter this effect, several finite-field methods have been 
proposed (e.g. Seitz & Schneider 1996). Secondly, since the 
shear arises from the gradient of the gravitational potential, 
the mass reconstruction is only known to within an additive 
constant. Accordingly, masses based on shear measurements 
are subject to a possible upward correction arising from a 
'mass sheet degeneracy'. 
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With sufficiently wide-field data it is possible to make 
the assumption that the surface mass density will approach 
zero at large distances from the cluster. However, with in- 
dependent knowledge of the magnification of the lens, the 
mass-sheet degeneracy can be broken regardless of the field 
of view of the data. Two methods have been proposed to 
make use of this magnification information and calibrate 
the absolute scale of the mass distribution, either through 
the change of image size at fixed surface brightness (Bartel- 
mann & Narayan 1995) or source counts (Broadhurst, Taylor 
& Peacock 1995). 

This paper is concerned with exploring the role that in- 
frared imaging offers in the gravitational depletion (or 'con- 
vergence') method for estimating the total masses of clus- 
ters. The depletion method was first suggested by Broad- 
hurst, Taylor & Peacock (1995) who predicted the diminu- 
tion in background galaxy surface number density as a func- 
tion of radius expected behind a lensing cluster. Here we are 
concerned with extending the original test to near-infrared 
wavelengths where, in principle, there are significant advan- 
tages, namely the fiatter number-count slope and a more ac- 
curate colour discrimination between foreground and back- 
ground populations. The unique wide-field capabilities of 
the panoramic near-infrared Cambridge Infrared Survey In- 
strument (CIRSI, Beckett et al. 1998) allow us to test the 
method on the rich cluster AbeU 2219 (z=0.22). 

A plan of the paper follows. In §2 we review the gravi- 
tational depletion method illustrating the difficulties associ- 
ated with its implementation at optical wavelengths and the 
potential gains of repeating the experiment at near-infrared 
wavelengths. In §3 we present new observations of AbeU 2219 
made at the prime focus of the 4.2m William Herschel tele- 
scope and discuss the techniques used to reduce the data as 
well as the methods used to create a sample of background 
galaxies. §4 discusses the depletion signal observed in the 
context of various mass models and reviews the uncertain- 
ties involved. In §5 we discuss the prospects of routinely es- 
timating cluster masses using this method both with CIRSI 
and with the upcoming suite of wide field infrared survey 
telescopes. 



2 GRAVITATIONAL DEPLETION 

The gravitational depletion or the 'convergence' method of 
breaking the mass-sheet degeneracy (Broadhurst, Taylor, & 
Peacock 1995) relies on the change in the surface number 
density of background galaxies induced by the magnification 
effect of a gravitational lens. Since only source counts are in- 
volved, exquisite image quality (essential for shear measure- 
ments) is not necessary. Furthermore, as the effect depends 
on the magnification /x, absolute mass estimates are possi- 
ble if the redshift distribution of the background sources is 
reasonably well-understood. 

The intrinsic (unlensed) counts no of galaxies brighter 
than some limiting magnitude m are transformed to the ob- 
served (lensed) counts n by 

n{< m) — no{< m)ii"^°'~^ , (1) 
where a is the logarithmic slope of the number counts, 

_ dlogno(m) 

dm ■ (2) 



Two competing effects serve to change the lensed surface 
number density of the background galaxies. Source mag- 
nification clearly increases the surface number density by 
magnifying galaxies that would otherwise be fainter than 
the limiting magnitude. However, focusing within the beam 
dilutes the overall surface number density. The net effect 
depends on the value of a: for a > 0.4 there will be an over- 
all increase in observed surface number density, while for 
Q < 0.4 a depletion is measured. 

Unfortunately, at the limits where a sizeable fraction of 
field galaxies are expected to be behind an intermediate red- 
shift cluster, the slope of the optical field counts, a « 0.4, 
produces only a weak effect. In order to demonstrate the 
effect in one of the most massive clusters known, AbeU 1689 
(2=0.18), Taylor et al. (1998) restrict their analysis to a red 
subsample known a priori from blank field studies to have 
a fiatter slope (a < 0.4). By colour-selecting sources red- 
der than the sequence of cluster spheroidals, Taylor et al. 
simultaneously secure a background population whose a is 
sufficiently low for the depletion method to work, and with 
an unlensed surface number density of 12 arcmin"^. Fort, 
Mellier & Dantel-Fort (1997) also search for an optical de- 
pletion effect behind the cluster CL0024-)-1654, by restrict- 
ing their search to the magnitude ranges 26 < B < 28 and 
24 < / < 26.5 where the slopes are found to be qb = 0.17 
and Qfi = 0.25. However, they apply no colour selection to 
remove faint cluster members in this magnitude range, and 
furthermore are left with an extremely low number density 
of only a few galaxies per square arcmin. Similarly, Athreya 
et al. (in preparation) use photometric redshifts to isolate 
a background population for their weak lensing analysis of 
MS 1008-1224, but are hampered by a small field of view 
and background clustering. Clearly, balancing accurate dis- 
crimination of the lensing foreground from the background 
population while maintaining both a fiat slope and sufficient 
numbers of background galaxies is of crucial importance, 
and it is for this reason that we turn to wide-field infrared 
observations. 

Here we are concerned with extending the depletion 
method to near-infrared wavelengths, and we illustrate the 
possible advantages via an initial application to the rich 
cluster AbeU 2219. The slope of the number counts fiat- 
tens at longer wavelengths because of a reduced sensitivity 
to intermediate-redshift star forming galaxies (Ellis 1997). 
Moreover, this flattening occurs at progressively brighter ap- 
parent magnitudes. Without any colour selection, the slope 
of the counts at infrared wavelengths is sub-critical with 
a ^ 0.25 at a relatively bright magnitude of A' ~ 18. This 
contrasts with optical counts which flatten only at very faint 
limits (e.g a ~ 0.40 for V = 22 - 24 flattening to a = 0.28 
beyond V = 24 (Small et al. 1995b)). Moreover, red-infrared 
colours (e.g. I — K) are more sensitive to redshift than to 
spectral class in the redshift range of interest. The opposite 
is true for say, B—I (see Fig. |^) . This leads to a much cleaner 
and efficient method of eliminating likely cluster members: 
we retain both the flat slope and sufficient numbers required 
for measurement of the depletion effect. The major draw- 
back of measuring depletion signals in the near-infrared un- 
til now has been the absence of panoramic infrared detec- 
tors capable of surveying large areas of sky rapidly to the 
required depth. With the commissioning of the Cambridge 
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Infrared Survey Instrument (CIRSI, Beckett et al. 1998), 
this becomes a practicality. 



3 DATA 

3.1 Strategy 

A number of factors enter when considering the merits of 
undertaking depletion studies at near-infrared wavelengths. 
Foremost, we can expect the counts to flatten considerably 
at cosmologically-significant depths for any passband long- 
ward of 1/im. Whereas classical number-count studies have 
been almost exclusively undertaken in the K- or Jf'-band 
(e.g. Moustakas et al. 1997 obtain a = 0.23 for 18 < A" < 23; 
Gardner et al. 1993 obtain ct = 0.26 for K > 18), recent 
counts at H (Yan et al. 1998) show httle change in slope, 
with Q = 0.31 ± 0.02 for 2Q < H < 24.5. Our lensing study 
will be undertaken using the //-band filter. 

Secondly, in terms of colour-selection, the degeneracy 
between redshift and spectral class likewise improves dra- 
matically when infrared magnitudes are added, particularly 
when account is taken of likely photometric errors. Fig. ^ il- 
lustrates a typical measurement of the I—H and B — I colour 
of a cluster early- type sequence at intermediate redshift. For 
the optical-optical colour we see that the single measurement 
is not enough to break the degeneracy between colour and 
redshift: an object bluer than the cluster sequence may be a 
low redshift elliptical or a late-type galaxy at any redshift. 
But when the optical-infrared colour is considered, the sensi- 
tivity to spectral type is greatly reduced, and we see that the 
bluer and redder galaxies map more cleanly onto foreground 
and background populations. One can therefore select more 
red objects with redshifts greater than that of the cluster. 
This accurate and efficient discrimination between the two 
populations is of great importance for our depletion study, 
which depends both on having a flat number count slope and 
a sufficient surface number density of background galaxies. 

Our principal goal is to explore the depletion expected 
in Abell 2219 in the context of earlier studies of this cluster. 
This is an X-ray luminous (Allen et al. 1992) cluster lying at 
a redshift of 0.22, with an X-ray temperature of Tx ~ 11.8 
keV and showing no evidence of a cooling flow (Allen & 
Fabian 1998). The distribution of the X-ray gas is elliptical 
and misaligned with the cD galaxy on small scales. Small 
et al. (1995a) report two systems of giant arcs with unde- 
termined redshifts: a 'red' arc (possible two merging images 
of a single background source) and a 'blue' arc consisting of 
three separate segments. The red arc is shorter and brighter, 
and is prominent in our //-band image (Fig. |^ . 

3.2 Observations 

The Cambridge Infrared Survey Instrument (CIRSI, Beck- 
ett et al. 1998) is a wide-fleld infrared imager consisting of 
a mosaic of four IK x IK Rockwell Hawaii HgCdTe detec- 
tors. Each detector is separated by a gap of approximately 
one detector width so that multiple pointings can be easily 
arranged to produce a larger contiguously-imaged field. 

Table 1 summarizes the observations that form the 
basis of this paper. The //-band CIRSI observations of 
Abell 2219 were made at the prime focus of the 4.2-m 




0.5 1 1.5 2 



Figure 1. Dependence of optical- infrared and optical-optical 
colour on redshift for early- and late-type galaxies derived from 
the GISSEL96 spectrophotometric codes (Bruzual & Chariot, in 
preparation) with and without evolution. The shaded area shows 
a range of colour when taking into account typical photometric 
errors of ±0.2 for the measurement of a typical early- type clus- 
ter sequence at intermediate redshift. Clearly the addition of the 
infrared magnitude reduces the dependence of colour on spectral 
type and allows for a more efficient selection of a background pop- 
ulation. This is an obvious advantage over using optical colours 
alone to distinguish between foreground and background popula- 
tions. 

William Herschel Telescope (WHT); this offers a plate scale 
of 0.32 arcsec pixel"^. Each CIRSI detector has a field of 
view of 5.5 x 5.5 arcmin so the mosaic of four detectors gives 
an instantaneous total field of view of 11 x 11 arcmin. 

The cluster was observed on the nights 1 to 4 July 
1999 using a field configuration chosen to optimally overlap 
with existing B- and /-band EEV CCD data of plate scale 
0.24 arcsec pixel" ^. Two adjacent CIRSI pointings were ar- 
ranged to ensure two contiguous strips of approximately 
5 X 20 arcmin each, one of which overlaps the two-pointing 
16 X 16 arcmin EEV field (Fig. |). 

The optical CCD data were taken as part of a com- 
panion programme and its reduction is discussed elsewhere 
(Bezecourt et al., 2000). The CIRSI observations were taken 
using short (30 sec) exposures with a 9 point pattern 
dithered with an offset of 15 arcsec. Larger offsets of the 
order of 50 arcsec were made before starting a new dither 
sequence. Typically 10 exposures were taken at each posi- 
tion, and the first exposure in each sequence was processed 
separately due to instability in the large scale bias level. 

Data reduction was done using routines from the CIRSI 
data reduction pipeline CIRDR (Chan et al. 1999, Hoenig et 
al. 1999) in iraf.^ For each dither sequence, a median sky 
image was created. The sky image was then scaled to and 
subtracted from each of the individual science images. The 
objects were located and a mask for each frame was con- 

t The Image Reduction and Analysis Facility (iraf) is dis- 
tributed by National Optical Astronomy Observatories, operated 
by the Association of Universities for Research in Astronomy, Inc., 
under contract to the National Science Foundation. 
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Table 1. Summary of observations 

Pointing a (J2000) 5 (J2000) total exposure time seeing 



CIRSI //-band, cluster centre 
CIRSI H-band, offset fields 
/-band EEV 



16 40 20.5 46 42 44.f 
16 39 50.2 46 42 44.! 
16 40 20.5 46 42 44.! 



5.3 h 
5.5 h 
1.0 h 



0.9" 
0.9" 
0.8" 



5 arctnin 



chipC: chip B: rhip R: 

easterncffset centralpoirting western o ff set 




EEV data 

Figure 2. The field of view of CIRSI after two adjacent pointings. The resulting two strips each measure ~ 20 X 5 arcmin. The dashed 
lines indicate the region which overlapped with the optical data. Treating the individual CIRSI chips as separate images, the western 
offset field, central cluster field, and eastern offset field (all with overlapping optical data) are labelled. North is up, east is to the left. 



structed using routines from Pat Hall's phiirs infrared data 
reduction package (Hall, Green & Cohen 1998). The origi- 
nal images were medianed a second time using the object 
masks to prevent contamination, and this second-pass sky 
was again subtracted from each of the science images. 

The resulting sky-subtracted images were averaged us- 
ing integer pixel offsets and a 3-sigma clipping rejection al- 
gorithm was applied. To construct the final image, the com- 
bined dither sequences from each night's observing were av- 
eraged in the same manner. Finally, the noisier regions at the 
edge of each mosaic where the overlap of dither sequences 
was incomplete were trimmed. Treating each detector as a 
separate image, one dither sequence was calibrated using 
measurements of a standard star from the UKIRT NIR cat- 
alogue (Hawarden et al. , in preparation). A set of secondary 
standards within the science image were then measured and 
used to calibrate the final mosaiced image, taking care to 



consider possible discrepancies in the overlap regions be- 
tween the chips. 

We used SExtractor 2.0 (Berlin & Arnouts 1996) to de- 
tect sources and perform isophotal photometry, selecting a 
detection and analysis threshold of 1.5ct above the sky back- 
ground. The //-band images were registered to the /-band 
data and SExtractor was run in double image mode. This 
allowed us to detect the objects in the /-band and measure 
both magnitudes with the same isophote. The catalogues 
were visually inspected and spurious sources removed. 

Fig. ^ shows the final reduced //-band image and in- 
dicates the degree of overlap with the EEV CCD data of 
Bezecourt et al. (2000). 
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Figure 3. Comparison of tlie CIRSI /f-band number counts witli 
those from tlie literature. Shown are the two CIRSI offset fields 
(covering 27 arcmin^ each), the Yan et al. (1998) //-band survey, 
and the Gardner et al. (1993) and Moustakas et al. (1997) K- 
band surveys. For comparison purposes the K-hand points have 
been shifted by one magnitude according to the approximation 
H — K = 1. The counts from the literature have all been corrected 
for incompleteness at faint magnitudes. 



3.3 Completeness and Contamination 

To understand the magnitude limit of our //-band data, 
we performed simulations by creating noise-only images: we 
randomized the offsets and thus created a misaligned mosaic 
with the same dimensions and noise properties as the science 
images, but with all the objects removed. We added artifi- 
cial galaxies of known magnitudes to these images using the 
artdata package in IRAF and recovered them using SEx- 
tractor. We found that for a detection in the //-band, the 
50% completeness level (the magnitude at which Sextrac- 
tor recovered 50% of the input galaxies) agreed well with 
the turnover in the observed number counts at // = 21.7. 
However, for purposes of measuring colour, we made the 
detections in the /-band and so were able to push the in- 
frared magnitude limit even deeper. Fig. ^ shows our result- 
ing //-band number counts compared to previous infrared 
surveys. 

In addition to adopting a strictly IR-limited sample, 
we have also explored alternative detection strategies not- 
ing that we seek the relative depletion in the counts as a 
function of position, i.e. uniformity of detection is more im- 
portant than completeness. We show in Appendix A that, 
for the purposes of our depletion measurements, it is not 
necessary to have a complete catalogue, so long as the in- 
completeness function responsible for the observed fall-off 
from the intrinsic power-law distribution is the same for the 
offset fields as for the lensed fields. As we determined the 
offset counts simultaneously under identical observing con- 
ditions, this holds true for our sample. We therefore push 
the limits of the //-band catalogue past the completeness 
threshold and adopt a working limit of // = 24 for our sam- 
ple, so that it remains //-limited. 

As discussed in Appendix A, it is important to con- 
sider not only the completeness of the sample but also its 
contamination by spurious sources arising from noise peaks. 



23 < /f < 24 

24 < // <30 
H =INDEF 



18 
235 
247 



3.6% 
47.0% 
49.4% 



Table 2. Noise fraction for the //-band image. Magnitudes were 
measured in 500 randomly placed apertures on a pure noise H- 
band image. The results show that the number of false objects 
brighter than our magnitude limit of // = 24 is negligible. 




H mag 

Figure 4. Colour-magnitude diagram for all objects in the clus- 
ter field. The cluster galaxies are indicated by filled circles; the 
population of objects selected for lensing purposes are those with 
colours redder than the cluster sequence and with 18 < // < 24 
(dashed line). 



especially when the counts are extended past the complete- 
ness limit. If we define the noise fraction as the fraction of 
false objects per magnitude bin, we may proceed with an 
incomplete sample only if the noise fraction brighter than 
our limiting magnitude is small. To estimate this we used 
the same pure noise image discussed above, and measured 
the magnitudes in 500 randomly placed apertures. Table 2 
shows the resulting distribution of magnitudes, with only 
3.6% of 'detections' falling below the H — 24 cutoff. This 
should be considered as an upper limit, since the detection in 
our sample is performed in the deeper /-band image. More- 
over, contamination, if it is left unaccounted for, will tend to 
reduce the lensing signal. We therefore conservatively ignore 
this small effect in our analysis. 



3.4 Removal of Cluster Members 

An essential precursor in examining a possible depletion ef- 
fect in Abell 2219 is the location of the lensing and unlensed 
population of galaxies, the latter of which is, of course, dom- 
inated for such a low redshift system by cluster members. 
The sequence of cluster galaxies is clearly visible on the I—H 
colour magnitude diagram shown in Fig. ^. Taking into ac- 
count known members and photometric errors, cluster galax- 
ies can be optimally excised according to the relation 

\I-H- (3.266 - 0.0865//)! < 0.2. 
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H mag 

Figure 5. Number-magnitude relation for non-cluster galaxies in 
the offset fields. Open circles are galaxies redder than the cluster 
sequence; filled circles are objects bluer than the cluster sequence. 
Note that the slope of the counts for the blue population is not 
significantly steeper than that of the red population (as is seen 
in the optical, e.g. Taylor et al. (1998)). Nevertheless we consider 
only the red objects for lensing purposes in order to minimize 
contamination by foreground galaxies. 



We can divide the remaining population of field galaxies 
into 'red' and 'blue' (unlensed) field populations, adding the 
further constraint that 18 < // < 24 and noting that the 
blue population could be composed of both foreground and 
blue cluster members. 

The two CIRSI chips on either side of the central cluster 
pointing provide us with flanking fields to characterize the 
field population. Combining the catalogues for the east and 
west flanking fields, we apply the same colour-selection to 
determine the properties of the unlensed background pop- 
ulation. Fig. ^ shows the resulting number counts for the 
sub-populations, which are described by 

log Nred = (0.185 ± 0.017)1/ + (0.23 ± 0.35), 

and 

log Nuue = (0.222 ± 0.030)iy + (-0.89 ± 0.61). 

for the range 19 < < 22 that is well described by a 
power-law distribution. Up to H = 24, the number den- 
sities of the background populations are no — 20.3 and 
no = 10.5 arcmin"^ for the red and blue populations, re- 
spectively. 

Finally, to account for obscuration of background galax- 
ies by foreground cluster members, we created a mask im- 
age using the SExtractor 'segmentation' option. This assigns 
each pixel in the image a value according to the object that 
contains it. We then used the location of the cluster galaxies 
(as well as those few non-cluster objects with H < 18) se- 
lected above to extract those pixels and form a mask which 
we used in the subsequent depletion analysis. 



4 RESULTS 

Our plan will be first to demonstrate the existence of the IR 
depletion behind Abell 2219 using the I — H colour diagram 



blue arc 




Figure 6. The central 2.7 X 2.7 arcmin of the _ff-band CIRSI 
image, showing the red arc reported in Small et al. (1995a). The 
second, bluer arc system is only barely visible in this infrared 
image. 



to isolat e cl uster members and select a background popu- 
lation (§3.4). We will then attempt to quantify the effect 
by employing maximum likelihood methods to fit simple, 
one-parameter models to the mass distribution for the clus- 
ter. Finally, we will consider the effect of uncertainties in 
the adopted mean surface number density on the maximum 
likelihood results. 



4.1 Depletion Analysis 

Fig. 1^ shows the radial density profile for field galaxies se- 
lected as discussed in §^^, taking into account the area 
masked by the colour-selected cluster galaxies (3.0 arcmin^ 
in total out of the 27.5 arcmin^ central field). We uti- 
lize both fianking fields in order to determine the un- 
lensed surface number density of the red population to be 
no=20.3 arcmin"'^. With no colour selection this density 
rises to nT=43.2 arcmin"^ which gives some indication of 
the likelihood that non-cluster galaxies might be discarded 
via this colour cut. 

Fig. gives a clear detection of the depletion signal at 
the 3(T level within a diameter of ~80 arcsec (~350 kpc for 
h — 0.5, Qo = 0.5). However, its absolute significance de- 
pends critically on the adopted value of no. We estimate 
the uncertainty in our measurement of no (indicated by the 
dashed lines in Fig. |^ using the dispersion in counts in the 
full //-band dataset shown in Fig. ^. The effects of this un- 
certainty are discussed further in §4.3. 

To investigate the magnitude and significance of this 
possible depletion, we adopted a maximum likelihood ap- 
proach based on that developed by Schneider, King & Er- 
ben (2000; hereafter SKE). In this formulation, which avoids 
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r (arcmin) 6^ (arcscc) 



Figure 7. Ratio of observed number density to unlensed number 
density of red galaxies with distance from the cluster centre, using 
a bin width of 14 arcsec. All measurements take into account the 
area in each bin obscured by the cluster galaxies. A depletion of 
number counts is clearly seen at small radii. The dashed lines 
indicate the errors in the measured value of no, the unlensed 
number derisitv (the effect of this uncertainty is considered in 
Section j }4.3| ). The thick line shows the expected depletion profile 
induced by a SIS with (t„ ~ 800 km s~^. 



Figure 8. Maximum likelihood analysis for depletion of the 
red background galaxy sample (solid circles). The model corre- 
sponding to the peak is a SIS parametrized by the Einstein ra- 
dius 9^ = 13.7 arcsec. By contrast, the same analysis for the red 
galaxy sample in one of the offset fields (open circles) peaks at 
9^ = arcsec and shows no evidence of depletion, as does the 
sample of blue (foreground) galaxies in the cluster field (open 
triangles). The curves have been vertically shifted to the same 
zeropoint for clarity. 



the loss of information induced by the radial binning (as in 
Fig. the log- likelihood equation takes the form: 



+ (2.5a - 1 



N 

E 



ln^(6'i 



(3) 



where no is the unlensed number density of background 
galaxies (arcmin"^), fi is the magnification, 9i is the posi- 
tion vector of the ith galaxy in the field with respect to the 
cluster centre, A'' is the total number of galaxies observed, 
and Q = dlog N/dm is the intrinsic logarithmic slope of the 
number counts. The first term in the log-likelihood function 
addresses the probability of finding A'^ galaxies in the field of 
view given the lens model /i(0) and the population parame- 
ters a and no , while the second term concerns the probability 
of finding each galaxy i at position 6i. By maximizing I we 
find the most likely parameter(s) for a given lens model. In 
the subsequent sections, we examine two single-parameter 
mass models. 

The derivation of this expression is presented in Ap- 
pendix A. Importantly, this appendix also generalizes the 
likelihood method of SKE to include the efi'ect of incom- 
pleteness. We showed that, as long as the intrinsic unlensed 
counts follow a power law, incompleteness is very simple 
to account for: one simply uses the observed (incomplete) 
unlensed density no in place of the intrinsic (complete) un- 
lensed density no in the likelihood function. This allows us 
to include, in our sample, galaxies which are fainter than the 
completeness limit and therefore of improving our signal-to- 
noise, without fear of introducing a bias. 

4-1.1 Model 1: Singular isothermal sphere 

First, we model the cluster according to a singular isother- 
mal sphere (SIS) parametrized with the Einstein radius, 9e: 



Msis(S) = 



1 - es/e 



(4) 



The filled circles in Fig. g show the resulting log- 
likelihood curve. The sharp dips are a result of the contribu- 
tion of the second term of Equation ^ (concerning the galaxy 
positions) to the log-likelihood function. As the likelihood is 
calculated for increasing values of the input 8e, a galaxy 
may happen to lie on the critical line {6i = Be). Since the 
magnification (and hence the depletion) is formally infinite 
at this radius, the probability of finding a galaxy vanishes, 
and the method rejects this particular model. This results 
from the fact that we have effectively modelled the galaxies 
as point sources. In practice, galaxies are extended objects 
and therefore are not subject to infinite magnifications, but 
instead are stretched into giant arcs on the critical line. In 
our analysis, we simply identified the peak as the absolute 
maximum, without any interpolation. As we will see below, 
this does not introduce any observable bias in the determi- 
nation of the model parameters. 

The peak of the likelihood function is reached at 
9e = 13.7 arcsec, which is consistent with the location of 
the red giant arc located 13 arcsec from the cluster centre 
(shown in Fig. 

We perform a test of the depletion by applying the same 
maximum-likelihood test on one of the offset fields. Using 
the same colour selection criteria to define a population of 
red objects we find no evidence for any depletion efi'ect. Fur- 
thermore, we also test the population of 'foreground' galax- 
ies (those with colours bluer than the cluster sequence) in 
the cluster field, and again find as expected no evidence for 
depletion. The log-likelihood functions for these two sam- 
ples, both peaking a.t 9e ~ arcsec, are also shown as open 
symbols in Fig. 
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4.1.2 Model 2: The NFW profile 

An alternative model to the SIS is the universal density 
profile for dark matter halos proposed by Navarro, Frenk 
& White (1997, 1996, 1995). Given the uncertainties in the 
data it is unlikely that we could use this technique to distin- 
guish between two models. It is nevertheless of interest to 
demonstrate how the method could be applied to a second 
model. 

The NFW density profile follows 

ScPc 



p{r) = 



(5) 



(r/r,)(l+r-/r,)2' 

where pc is the critical density at the redshift of the lens. 
The two parameters of the model are contained in the scale 
radius and the concentration parameter c, so that the 
characteristic overdensity is 



200 



3 ln(l + c)-c/(l + c). 



(6) 



In contrast to the SIS, this profile fiattens towards the 
core and is capable of producing radial as well as tangential 
arcs (Bartelmann 1996). For a projected radial distance R, 
we can define a dimensionless distance, x — R/vs. Wright & 
Brainerd (1999) give the formulation for the radial depen- 
dence of the surface density E(x) and the shear 7(2;). The 
convergence is then simply n{x) = E(a::)/Ecrit, where the 
critical surface number mass density 

" 4-KGDdDds ^ ' 

depends on the angular diameter distances Ds,Dd, Dds from 
observer to source, observer to lens, and source to lens, re- 
spectively, and the speed of light v^. The magnification is 
then simply 



MNFw(a;) = 



(l-K(a;))2-7(a;)2) 



(8) 



(e.g. Schneider, Ehlers and Falco 1992). 

Using /iNFW in Equation we perform the maximum 
Ukelihood analysis for an NFW profile on Abell 2219. We 
reduce the model to one parameter by assuming a reasonable 
value for the concentration parameter, namely c = 6, and 
fitting for Vs. We find a maximum at Vs corresponding to an 
angular length of 52 arcsec. At the redshift of the cluster this 
gives a scale length of 125/i~^ kpc (where Ho = 100^ km s~^ 
Mpc~^), which compares well to the ~ 250/i~^ predicted 
from numerical simulations of CDM models (Bartelmann 
1996). 

Fig. ^ shows the form of an NFW depletion curve result- 
ing from the substitution of /xnpw in Equation |l[ for such 
a cluster with c = 6 and Vs = 125^"^ kpc . The observed 
depletion curve of Abell 2219 is overplotted for comparison, 
and we see that the data are also consistent with this model. 
However, the data do not allow us to distinguish between the 
NFW and SIS models. In the further analysis we choose to 
continue with the SIS for simplicity. 

4.2 Simulations 

To further test the robustness of the depletion effect we have 
simulated a population of galaxies with the same proper- 
ties as our putative background population. In our observed 




r (arcmin) 

Figure 9. The radial depletion curve for a NFW cluster with 
Ts = 12bh~^ kpc and c = 6. The inner dip corresponds to the 
radial arc at Xr = 2 arcsec while the outer dip corresponds to 
the tangential arc at = 5 arcsec, for a source population with 
a = 0.185. Overjayed are the points the Abell 2219 depletion 
■ Fig. I 



curve, as in 



sample the intrinsic power law form of the background num- 
ber counts is modified by the incompleteness function. As 
shown in Appendix A, it is however equivalent to use a renor- 
malized power law for the counts in the simulations. Specifi- 
cally, we keep the same observed slope, a = 0.185 but lower 
the normalization to maintain the same observed unlensed 
number density fio = 20.3 arcmin"^. 

Scattering the simulated galaxies randomly across the 
field of view, we then adjusted their positions and magni- 
tudes as if they were lensed by a singular isothermal sphere 
with = 13.7 arcsec. We created a catalogue of those galax- 
ies whose lensed magnitudes were brighter than our limiting 
magnitude and whose lensed positions were not obscured 
by the Abell 2219 cluster mask and applied the same maxi- 
mum likelihood technique. Fig. ^ shows the distribution of 
the resulting best-fitting models for 200 realizations. We see 
that the distribution of recovered Einstein radii is centred 
around the input value of 13.7 arcsec, and deduce the 95% 
confidence interval to span the range 6.8 < 9e < 26.5 arcsec. 

To translate the estimates of into an statement re- 
garding the cluster mass, we use the relation for a SIS: 



(9) 



where is the velocity dispersion of the lens, Dds and Ds 
are the angular-diameter distances from lens to source and 
from observer to source, respectively. 

Modelling the background N{z) from the K-hand sur- 
veys of Cowie et al. (1996), we find the median redshift 
for the background is (z) — 1.0 for H < 22, increasing to 
(z) — 1.3 at _ff < 24. Given zi^ns = 0.22 and assuming a cos- 
mology of {Ho = 50 km s ~^ Mpc go ~ 0.5} we obtain 
cr„ = 842t\ll km s"^ for (z) = 1 and cr„ = SUtlH km 
for (z) = 1.3. The quoted errors refer to the 66% confidence 
limits derived from the above simulations. Clearly the result 
is more sensitive to the uncertainties in the measurement of 
6e than the uncertainties in the background redshift distri- 
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^ with cluster mask 
Q no cluster mask 



9, (arcsoc) 



Figure 10. Distribution of the recovered Se from 200 realiza- 
tions of simulated lensing by an isothermal sphere, using the ob- 
served background population parameters (no = 20.3 arcmin"^, 
a = 0.185). The darker shaded region indicates the resulting 66% 
confidence interval in the range 11.1 < 6e < 17.06 arcsec, cen- 
tred around the input value of 8^ = 13.7 arcsec. The lighter 
shaded region indicates the 95% confidence interval in the range 
6.8 <9e < 26.5 arcsec. 



bution. Therefore, we consider the background to lie on a 
single sheet at the median redshift, and placing the galaxies 
at 2 = 1.3 and adding the errors in quadrature we obtain our 
final result of (t„ = SM^^Jg km s~^ (66% confidence level). 
Our SIS model is in good agreement with the result from 
Small et al. (1995a), who find cr„ ~ 930 km s~^ using strong 
lensing constraints. Note that the precision of the strong 
lensing constraint would be tighted if the redshift of the red 
arc were known. Our results are also consistent with the 
optical weak shear analysis of Bezecourt et al. (2000), who 
find (Ti, = 1075 ± 100 km s~^ from an optical weak shear 
analysis. In principle, by adding one additional parameter 
to our model we could also constrain the amplitude of a 
mass-sheet, and therefore break the degeneracy inherent in 
reconstructions involving shear only. However, fitting an ex- 
tra parameter is not warranted by the data, especially after 
taking into account the uncertainty in no (see §4.3 below). 

The above simulations only enable us to deduce the er- 
ror on 9e assuming our initial estimate of the model param- 
eters is correct. To test the method across a wider range of 
input lenses, we have simulated lenses with input parameters 
6e ~ 0, 10, 20, and 30 arcsec, and performed the maximum 
likelihood analysis 25 times in each case. Fig. |l^ shows the 
recovered Oe vs the input 9e, both with and without the 
cluster mask. We see that neither the dips in the likelihood 
function, nor the masked area introduce any significant bias 
into the method. In addition, we find that obtaining a best- 
fitting value of 13.7 in the absence of lensing, i.e. a false 
detection, is unlikely, and we reject this null hypothesis at 
the 6cr level. This figure provides a measure of the signifi- 
cance of our detection of the lensing signal if the uncertainty 
in 710 is ignored and the assumed model is correct. 



Figure 11. Simulations showing recovered Be vs input 6e, with 
and without a cluster mask (the open circles are offset slightly in 
the +x direction for clarity). The cluster mask, which obscures 
about 10% of the field including much of the central region, has 
no significant effect on the recovered 9e (aside from increasing 
the errors), even given the low numbers of background galaxies. 



4.3 Effect of uncertainty in number counts 

SKE demonstrate that, to first order, most of the magni- 
fication information is provided by the unlensed number 
density no. Without adequate knowledge of this normal- 
ization, the mass-sheet degeneracy cannot be broken. Prior 
knowledge of the uncertainty in no can be included in the 
maximum-likelihood analysis. Assuming that the true no fol- 
lows a Gaussian distribution with mean fio and dispersion 
CTn = r?no then, as shown in SKE, the log-likelihood function 
becomes: 



= -no74-(/3-l)^ln^(6'0+A^ln 



no' 



(no — UoY 
2 {rjUof 



,(10) 



where / = / £'e[^i[6)]'^'^ and (5 = 2.5Qf. This can be maxi- 
mized with respect to no, yielding 

1^ = i (1 - r^'nj) + y i (1 - n^nj) + rj^N. (11) 

Substituting this value for no into Equation ^ yields the 
new log-likelihood function. 

To determine the value of rj for our own sainple, we 
turn to the entire H-hand dataset pictured in Fig. bl While 
the sample used in the depletion analysis is a colour-selected 
sample that incorporates the /-band data, the area in com- 
mon with the optical image allowed for only the two offset 
fields discussed previously. However, we attempt to estimate 
the uncertainty in no by sampling the counts in the entire 
//-band region. 

We divided the two //-band strips into four single- 
chip images each, discarding the chip containing the clus- 
ter and the two images taken with chip 3 (which was con- 
taminated by many hot pixels, leading to a high number 
of false detections). This left five offset images of equal 
area. We constructed //-limited catalogues on the central 
4.3 X 4.3 arcmin region of each image using SExtractor and 
examined the dispersion in the total //-band number counts 
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field N(18 <H <2A) 



468 
516 
483 
476 
394 



Table 3. //-band counts in 18.2 arcmin^ regions in five of the 
seven offset fields, yielding a fractional uncertainty of 8.6%. 



in the range 18 < // < 24 for each chip. The resulting distri- 
bution of raw counts, listed in Table 3, has mean N = 467 
and standard deviation a = 40, hence the fractional error 
is ?7h = 0.086. If the fraction of galaxies discarded during 
the colour selection to form the 'red' sample of background 
galaxies is constant from field to field (from the two multi- 
colour offset fields it appears to be ~ 50%), then this error 
will also apply to the subpopulation. 

We note that this is higher than the Poisson error for 
these number counts. However, if we consider the angular 
correlation function lo{9) of K-hand galaxies to A" = 21.5 
presented in Carlberg et al. (1997) and extrapolate to 
6 = 2.5 arcmin then cli{0) ~ 0.01. This yields a rough esti- 
mate of SN/N ~ ^/u;{2.5') ~ 0.1 which is in agreement with 
our measurement of the dispersion. Our uncertainty in the 
unlensed density is therefore dominated by clustering. 

We again turn to simula tion s to derive the error bars 
on our measurement. As in §4.2 we simulate a population 



that is lensed by a SIS with 6e = 13.7 arcsec, but in this 
case we draw the true value of no from a Gaussian distri- 
bution with mean 20.3 arcmin^ and dispersion cr„ = ryno, 
where rj = 0.086 as derived above. Using the log-likelihood 
function in Equation ^ and a fixed estimate no = 20.3, we 
show the resultin g d istribution of recovered 9e for 200 re- 
alizations in Fig, We see, by comparison with Fig. |l^ 
(the case where no was known precisely), that incorporating 
an error of ?7 = 0.086 produces a much broader distribution, 
with the result of greatly increasing the error bars on our 
measurement. Applying Equation ud to the real data, we 
find 9e = 1.4ti;4 arcsec (66% confidence) or Oe = lA^l^f 
(95% confidence). While the estimate of 9e itself is greatly 
reduced, the error bars make it compatible with the result 
of § [4. 2| within a 2a range. 

This section illustrates the principal weakness of the 
depletion method as pointed out in SKE, namely the vital 
importance of accurately measuring no. With a 5N/N ~ 0.1 
which is not unreasonable given clustering on these scales, 
the depletion effect is made more difficult to distinguish from 
variations in the background counts that are not due to lens- 
ing. However, this effect could be countered by selecting sim- 
ilar clusters (e.g. by their X-ray temperatures) and stacking 
the depletion signal accordingly, thus calibrating the cluster 
Tx-mass relation and obtaining an average cluster mass pro- 
file. Choosing a sample of 10 clusters would simultaneously 
increase the effective background surface number density by 
a factor of 10 while reducing the fractional error in no by a 
factor of \/To. Simulations show that for 10 clusters similar 
to Abell 2219, the 95% confidence region on 9e ~ 13.7 arcsec 
would then shrink from 0.9 < 9e < 29.6 arcsec (as quoted 
above) to only 7.7 < 9e < 16.1 arcsec; a feasible project for 
upcoming IR survey telescopes. Note that in practice scatter 
in the cluster properties would increase this range somewhat. 




Figure 12. The effect of uncertainty in the background num- 
ber counts from simulations. Shown is the distribution of the re- 
covered 9e for 200 realizations, with the true no drawn from a 
Gaussian distribution with mean 20.3 arcmin"^ and fractional 
error rj = HJi%. The lensed populations were analysed as in 
Equations hd and nil using the estimated no=20.3 arcmin"^. 
The 66% and 95% confidence limits are indicated by the dark 
(8.9 < 6»j5 < 18.9) and fight (0.9 <8e < 29.6) shaded regions. A 
comparison with Fig. hd shows that incorporating a significant 
uncertainty in the analysis increases the error bars on the mea- 
surement oi 9e- 



5 CONCLUSIONS 

We present a study of the depletion effect around Abell 
2219, the first done in the infrared, using the newly- available 
panoramic IR camera (CIRSI). We show (see Appendix A) 
that the sample can be effectively exploited beyond the com- 
pleteness limit, as long as the lensed field and the back- 
ground field obey the same incompleteness functions and 
have minimal contamination by false objects. This allows 
us to detect a clear dip in the radial number density pro- 
file of background galaxies at small distances from the clus- 
ter centre. The optical-infrared colours enable us to select a 
population of red background galaxies with a fiat number- 
count slope in order to optimise the lensing signal and reduce 
foreground-background confusion. 

For a population of red background galaxies with an 
extremely fiat slope (a = 0.185), we employ maximum like- 
lihood methods and a SIS model to derive an estimate of the 
Einstein radius 9e = 13.714 2 arcsec (66% confidence limit 
when uncertainties is no are ignored) , resulting velocity dis- 
persion (Tt, = 814li[3g km s~^. These values are consistent 
with the location of the redder of the two giant arcs and the 
estimate a„ ~ 930 km s~^ of Smail et al. (1995a). 

We examine the uncertainty in the number counts, and 
derive a fractional error of 8.6% on the normalisation of 
the background number density (consistent with cluster- 
ing on these scales). When this error is incorporated into 
the maximum-likelihood analysis the error bars become too 
large to make a precise statement about the magnitude of 
the lensing (although our previous measurement is not ruled 
out). This demonstrates the crucial importance of the back- 
ground density no for an accurate depletion analysis as dis- 
cussed in Schneider, King & Erben (2000). 
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Finally, while we cannot at present use our current data 
to distinguish between alternative models for the cluster 
mass distribution, we demonstrate the application of the 
method to the data using a second mass model. We show 
the typical depletion curve induced by a cluster with an 
NFW density profile, and using this profile as a second esti- 
mate of the cluster mass we derive a scale length of 125/i~^ 
kpc for the cluster. 

The wide field IR coverage available using CIRSI is an 
ideal match to the current generation of wide-field optical 
imagers. The availability of multicolour data taken simulta- 
neously but at large radii from the cluster centre provides 
offset information about the field population and eliminates 
the need to normalise the depletion curve to the values at the 
edge of the chip, and the infrared sample provides a cleaner 
discriminant between foreground and background and a flat 
number-count slope. 

The depletion method is an elegant and relatively sim- 
ple way to derive mass estimates of clusters of galax- 
ies. Deeper IR surveys will provide more accurate knowl- 
edge of oj{0) on arcminute scales, allowing us to bet- 
ter understand the degree to which clustering affects 
the method. Future studies of clusters with wide-field 
optical-infrared data (e.g. with the V ISTA telescope, 
tittp : //www-star . qmw . ac .uk/~jpe/vista| ) covering a wide 
wavelength range could provide more accurately selected 
background populations via photometric redshifts and allow 
us to add another sample of independent mass profiles to be 
compared with those derived from velocity dispersions, X- 
ray measurements, and strong and weak lensing. The prob- 
lem posed by the uncertainty in the background number 
counts can be overcome by selecting similar clusters (e.g. by 
their X-ray temperatures) and stacking the depletion sig- 
nal accordingly to obtain an average cluster mass profile; a 
feasible project for future IR survey telescopes. 
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APPENDIX A: EFFECT OF 
INCOMPLETENESS 

The likelihood function for the lensing depletion has been 
derived by SKE, for the case of a complete sample of back- 
ground galaxies. Here we generalize their results to include 
the effect of incompleteness in the sample. 

Let the unlensed differential counts of the background 
galaxies be denoted by where no is the number of 

galaxies per unit solid angle, and S is the flux. The lensed 
counts are given by 
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dn{S) _ 1 dnoiji) 



dS 



dS 



(Al) 



where ^ is the magnification. The observed (lensed) counts, 
in the presence of incompleteness, can be written as 



dh{S) _ dn{S) 



<t>{S), 



(A2) 



dS dS 

where (t){S) is the completeness fraction for a given flux 5*. 
Typically, 0=1 for large fluxes, and 4> <1 for small fluxes. 
In the following, we will continue to use the subscript o to 
denote unlensed counts, and a tilde to denote the observed 
counts. In practice, it is convenient to consider the total 
counts down to a detection limit S 



n(> S) 



dS 



, dn{S') 
dS 



(A3) 



In our case (see Fig. ^) and for many applications, the 
intrinsic unlensed counts are well described by a power law 
of the form 
driQ 



cx 5" 



no(> S) oc 5" 



(A4) 



dS 

where (3 is the flux slope parameter, which is related to the 
magnitude slope parameter a = d(log nn] /dm (Equation ^ 
hy 13 — 2.5a. In this case. Equation (Al) takes the simple 
form 

dn j dno 
Is / ~dS 

In this case, the observed lensed counts are given by 
dh 



n{> S) /no(> S) 



0-1 



(A5) 



ds^'^'' -ds- 



(A6) 



As a result, the mean integrated lensed counts (A'^) ob- 
served in the fleld are given by 



(N) 



dS 



dn 



where ho = /g 



dS 
dS'4>- 



ho / d'e [fi{e)Y 



(A7) 



is the observed surface number 



)S ■f' dS 

density of unlensed galaxies. 

The probability of finding the i**^ galaxies at position 9i 
is Pi oc h{9i) = no[/i(Si)]'^~^. Thus, the likelihood of finding 
A'^ galaxies in our incomplete survey at positions 9i is given 

by 



C = P{N-{N))\{ 



1-9-1 



(A8) 



where P{N; (N)) stands for the Poisson distribution with 
mean {N}. The log likelihood function is thus 



i^-ho / d''e[ii{e,)f''^ + {/3- 



N 



(A9) 



where constant terms have been dropped. 

These two expressions are identical to the corresponding 
expressions in SKE, except that incomplete counts (N) and 
no replace the complete counts (N) and no, respectively. 
Therefore, the effect of incompleteness is very simple to take 
into account. (We will drop the tilde in the rest of the text). 
In fact, there is no reason to restrict the sample to galaxies 
brighter than the completeness limit. On the contrary, it is 
desirable to include faint galaxies to improve the signal-to- 
noise ratio. 



Note that this simplification results from the factorisa- 
tion of n in Equation (A6) and thus only holds if the in- 



trinsic unlensed counts follow a power law. Note also that 
the likelihood £. is independent of the completeness function 
4i{S). Thus, two samples with different (ji's but with the same 
observed galaxy density no would be statistically indistin- 
guishable, as far as the depletion analysis is concerned. This 
aUows us to use a power law count with a rescaled normaliza- 
tion (corresponding to a constant (l>{S)) in the simulations, 
and thus to avoid modeling the completeness function. 

At the very faint end, the sample will not only be incom- 
plete, but it will also be contaminated by spurious sources 
arising from noise peaks. These spurious sources wiU be uni- 
formly distributed on the chip and will not be affected by 
lensing. While this contamination could easily be introduced 
in the likelihood function, it is hardly necessary in practice. 
Indeed, incompleteness typically becomes important before 
contamination does. As a result, a moderate magnitude cut 
ensures that the contamination fraction in small. In addi- 
tion, contamination tends to 'wash out' the lensing signal. 
From the point of view of detecting the lensing signal, it is 
therefore conservative to ignore this effect. 

This paper has been produced using the Royal Astronomical 
Society/Blackwell Science style file. 
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